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Abstract 

The high order corrections to renormalon are considered. Each new type of 
. insertions into the renormalon chain of graphs generates the correction to the 

asymptotics of perturbation theory of the order of ~ 1. However, this series of 
corrections to the asymptotics is not the asymptotic one (i.e. the m-th correction 
does not grow like m\). The summation of these corrections for UV renormalon 
may change the asymptotics by factor N s . For the traditional IR renormalon the 
£NJ ' m-th correction diverges like (— 2) m . However, this divergency has no infrared 

origin and may be removed by proper redefinition of IR renormalon. On the other 
, hand for IR renormalons in hadronic event shapes one should naturally expect 

| these multi-loop contributions to decrease like (— 2)~ m . Some problems expected 

upon reaching the best accuracy of perturbative QCD are also discussed. 
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)£ ' 1 Introduction 

On 

The running coupling constant seems to be an inevitable companion of any nontrivial 
Ph. renormalizable field theory. In its turn the asymptotics of perturbation theory for any 

quantity calculated in the theory with running coupling in general is determined by 
renormalon |], |] . The renewed interest in this kind of asymptotic estimates have 
been demonstrated in last few years 0, EL [71 ^, [II], O, O. It results even in 



X 



attempts |L4| to use renormalon for calculation of experimentally measurable quantities. 
Finally, the recent explosion of activity in the renormalon business has been stimulated 
by the understanding that the former may be related with the anomalously large (~ 
Aqcd/Q) nonperturbative corrections in hadronic shape variables |15|, [16|, [18|, [17|, [19|, |20fl. 

The accurate determination of renormalon-type asymptotics turns out to be not so 
simple problem. It was recognised @, 10, 12, 22 1 that the overall normalisation of the 



renormalon could not be found without taking into account of all terms of the expansion 
of, e.g., the Gell-Mann-Low function. However, surprisingly up to now nobody have 
tried to sum up this series of corrections to the renormalon. 

Therefore the main aim of this paper will be to discuss the possible consequences of 
summation of contributions from the arbitrary high order insertions to the dressed gluon 
line. As we demonstrate both by diagrammatic consideration and by direct analytical 
calculation, each new type of insertions generates the correction to renormalon of the 
order of ~ 1. However, the fc-th correction to the asymptotics for large k is not expected 
to have any k\ enhancement. Thus at least the series of corrections to the amplitude of 
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renormalon asymptotics is not the asymptotic series. On the other hand, the summation 
of this series leads to the sufficient change of the high order behaviour of the usual 
perturbative series. We show that for the UV renormalon taking into account of the 
high order corrections may change the iV-th term of the perturbation theory by the 
factor N s (although with S probably unknown even for QED). The m-th correction to 
the IR renormalon diverges like (— 2) m . Due to this divergency the IR renormalon does 
not exist in the usual sense. We propose however the proper rearrangement of the series 
which allows to avoid this difficulty. The IR renormalon now will be associated with 
certain non-Borel-summable series but the coefficients of this series themselves are the 
Borel-summable series in a s . Finally, the event shape renormalons, which are the most 
powerful among other renormalons, may be quite stable under taking into account of 
the high order insertions. The corresponding m-th correction is expected to decrease 
like (-2)" m . 

In accordance with their name renormalons are connected with the running of the 
coupling constant. Therefore in considering of the high order corrections to renormalon 
one is faced with the problem of renormalization scheme independent definition of the 
coupling. This problem may be solved quite naturally in QED. That is why, in particular, 
all the diagrammatic examples of this paper will be of the QED-type. The generalization 
to the QCD case is usually done in some heuristic way like e.g. the so called 'naive 
nonabelization'. Nevertheless, the first attempts to find the renormalization scheme 
independent definition of the effective charge in QCD have been done recently |24|| . 
Also we will not usually make much difference between photon and gluon, thus making 
the tacit assumption that the difference between theories is hidden somewhere in the 
coefficients bo, fei, 62, • • ■ of the Gell-Mann-Low equation. 

The contribution of the diagrams with exchange of one soft gluon(photon) to some 
"physical" quantity which accounts for the infrared(IR) renormalon has the generic form 

Rm = I a(k)^ . (1) 

Jk^Q C/ 1 

The Feynman graphs corresponding to this quantity are shown in fig. 1. Similarly the 




a b 



Fig. 1 The renormalon graphs with exchange of one gluon. The internal gluon line will be dressed in 
the following figures. 
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exchange of very hard gluon(photon) generates the ultraviolet (UV) renormalon 

Ruv= I a(*)2j^ . (2) 
Jfc>Q k b 

Now only the diagram of fig. la contribute. During the last two years it was understood 
|12| , [22] that the traditional UV renormalon (0) with exchange of only one hard photon 
does not give the largest contribution to the asymptotics. The diagrams with exchange of 
at least two photons turns out to be much more important. We will return to discussion 
of this new UV renormalon later while now our consideration of the traditional UV 
renormalon is of mainly methodological importance. 

In (|T],0) we have written down the effective running coupling constant a(k) = a e ff(k), 
which (at least for QED, but see also J24J) is trivially connected with the transverse part 
of the gluon propagator. The function a(k) satisfies the RG equation: 

^ = b a 2 + ha 3 + b 2 a* + ... , x = In (Q 2 /k 2 ) . (3) 

It is to be noted here that we have fixed the renormalization scheme by considering the 
effective charge. Thus our coefficients b 2 , 63, ... are neither the free parameters, nor the 
known, e.g., for MS scheme, b 2 (MS), b 3 (MS) . At first stage one may neglect bi,b 2 , ... 



in 



Rm = - ra(x)e- 2 *dx= r 2dx=Y, ( N\ . (4) 

a Jo Jo 1 — b a x ^— p \ 2 J 

Here we have chosen some convenient overall normalization of the renormalon. We will 
consider now only the asymptotics of the perturbation theory and leave the issue of the 
nonperturbative ambiguity of the integral (ffl) due to the Landau pole to the very end 
of the paper. The integral ([]]) describes adequately the contribution of a certain chain 
of Feynman diagrams only for k <C Q . It is seen from (|J) that the main contribution 
to the iV-th order of the expansion comes from k 2 ~ Q 2 e~ N l 2 . Thus the renormalon 
contributions to the first few terms of perturbation theory are completely irrelevant. On 
the other hand, for sufficiently large Q a lot of terms of the expansion (H) come from 
the region A% CD <C k 2 Q 2 , where the effective charge is small and the perturbative 
approach for calculation of a e ff ([|) seems to be useful. 

Very similar calculation for the UV renormalon @ leads to the same result as (|J) up 
to trivial replacement ^ — ■> — bo- This means on the one hand that the UV renormalon 
asymptotics in general is much more important (in 2 N times) than the IR one. On the 
other hand one may see that depending on the sign of 60 the series associated with one 
of the two renormalons should be non-Borel summable. 



2 Diagrammatic examples 

The main part of the current interest in renormalons seems to be concentrated on the 
consideration of their relation with nonperturbative corrections to various observables. 
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b, Q 2 ln(Q 2 /k 2 ) 



b a ln(QVk z ) 




Fig. 2 The simplest chain of diagrams Fig. 3 The example of diagram with 
corresponding to the renormalization two - loop insertion into soft gluon line, 
of the soft gluon line. Each bubble 
generates the factor 6qq;o In (Q 2 /k 2 ). 



However, just this part of the renormalon business suffers mostly on the hard to es- 
tablish assumption about the enhancement of nonperturbative corrections compared to 
the typical uncertainty of renormalon resummation. The only free of phenomenological 
input use of renormalon is for the approximate (with 1/N as a small expansion param- 
eter, where N is the number of the term of the perturbative series) calculation of the 
high order contribution of the usual perturbation theory. Therefore, before passing to 
straightforward but rather formal manipulations with the RG equation let us illus- 
trate the role of complicated contributions to renormalon by the few explicit estimates 
of Feynman graphs. 

Moreover, the experience of this diagrammatic consideration will help us in the fol- 
lowing to distinguish, to what extent the factorial growth of the series is connected with 
the true infra-red(ultra-violet) physics or appears due to a simple combinatorics. 

In this section we will consider only the IR renormalon. Generalization for the UV 
one is straightforward. The fig. 2 shows the chain of diagrams corresponding to We 
show only the QED-type diagrams without gluon self-interaction. Each of the N bubbles 
from fig. 2 generates the factor &o«ohi {Q 2 /k 2 ) in the integrand of ([!]),(§)• The difference 
between QCD and QED may be thought to be hidden in the factor 6 0) accompanying 
the single bubble. 

Now let us replace two of the simple bubbles by the more complicated diagram of 
fig. 3. The two loop bubbles generate the factor b\a^ In {Q 2 /k 2 ) in the integrand, which 
has one power of large logarithm less (or one «o more) than the leading order contribution 
(H). However, a large combinatorial factor iV appears due to a number of permutations 
of the second order bubble among the simple bubbles, leading to 



Attain 



2\ - / / v>\ l N-2 /, \ N 



b a In ( Or 



k 2 

Thus we see that taking into account one second order insertion into the soft gluon line 
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leads to the correction of the order of one to the trivial asymptotics (J4|). 



Consider now the more complicated diagram of fig. 4 with dressing of the internal 
gluon line of the second order bubble. To this end it is natural to write down explicitly 
the last integration over internal momentum of the two loop diagram 



k' 2 



6o«o In 



'Q 2 



dq 2 



n + 1 



bia In 




6o«o In 




(6) 



Thus up to the overall factor ^j-j- the contribution of diagram of fig. 4 coincides with 



Fig. 4 The dressing of internal gluon line of the second order bubble by n simple bubbles. 

that of fig. 3. Summation over n naturally leads to ln(iV). Taking into account a number 

q 



Fig. 5 Three loop insertion with dressing of two internal gluon lines by the simple chains of bubbles. 
The summation over n\ and ni allows to compensate all extra as. 

of large bubbles of fig. 4 allows to exponentiate the correction 

(^) N i V !exp(|.n(i V )) = (^)"' J vti V !. (7) 

This is the generally recognised expression for the IR renormalon. Our argumentation 
up to this stage repeats the line of reasoning of the paper [^] . However, the argument of 
the exponent in (|7]) was found with the ~ 1/ ln(iV) accuracy and therefore the nontrivial 
overall factor as well as the function of N, weaker than iV 7 , may appear in (|7]). 

Now let us consider the three loop correction (fig. 5 with n\ — n<i = 0). This 
contribution generates the factor 62^0 In (Q 2 /k 2 ) in the integrand of (^). Thus here we 
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have two extra a which at first glance could not be compensated by one combinatorial 
N and hence the diagram of fig. 5 seems to generate only the ~ 1/N correction to 
renormalon. However, let us see, what happens if one dresses the internal gluon lines 
of the three loop diagram. Now summation over the number of trivial insertions ni,ri2 
gives: 

**«o In (%) £ *l N ~ n !~ n *: 2 ~ ha In (%) x {a,Nf . (8) 
V k J TH+U2+1 \k 2 ) 

Here the factor (n\ + n 2 + l)" 1 appears after integration over the internal momentum of 
the large bubble, while (N — n\ — n 2 — 2) accounts for the combinatorics. We see that 
after dressing of all gluon lines the three loop (~ b 2 ) diagram generates the correction 
to renormalon of the order of ~ 1. One can easily show that four loop (~ 63), five 
loop (~ 64) etc. diagrams generate the corrections of the same order of magnitude. 
Previously the analogous proof of the importance of the high loop corrections was done 
by Mueller pl| but this result was not published. 



3 Upon summation of the corrections for the UV 
renormalon 

It is easy to integrate formally the renormalization group equation (|3|) 

1 1 b\ ( a \ 2 2 
\ — In I — - c 2 [a - a ) - c 3 (a - a ) - ... = b x , (9) 



a ao 60 veto' 

where, c 2 = b 2 /b — b\jb\ , c 3 = {b 3 /b — 1b 2 b\lb\ + b\/bl)/2, ... . For arbitrary k one 
has Cfc = {bk/bo — ...)/{k — 1). This "exact" solution still is too informative for us. First 
of all, most of the terms containing a in (|), namely c 2 a , c 3 «q, c 4 q;q, ... will contribute 
only to the ~ 1/N corrections to the asymptotics. Therefore we may write 

a = ^0 ^ 1Q ^ 

1 — 60^0^ — (bi/bo)ao\n(a/ao) — c 2 aao — c 3 a 2 ao — ... 

Let us introduce now the truncated running coupling 

= ^0 

' 1 - b a x - (6 1 /6 )a ln(at/a ) 

Now one may again formally expand the exact a (0) in the series 

a = a t (l + I3 2 a 2 t + foal + ...) , (12) 



The explicit formulas for a first few coefficients fo may be easily found. However we 
will be interested only in the large-Zc asymptotic behaviour of fo- Naturally for our 
definition of the effective charge (|j) the coefficients bk themselves form the asymptotic 
series bk ~ k\. Due to that one easily finds the asymptotics of (3k in ([12]). 



A = &( 1 + °(s» • (13) 
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In the previous section we have shown diagrammatically that contributions to the 
asymptotics from the high order terms of the RG equation (|]) ~ b 2 ,bs,b4,... (or now 
~ /3 2 , P3, @4, ■■■ dUt )) are not small. However, the contribution induced by the second 
(two loop) term b% a 3 play an outstanding role due to the additional enhancement by 
ln(iV). This is the reason for taking the truncated coupling at (0) as the expansion 
parameter in (|12D . 



The expression ( |TT|) itself is the transcendental equation for function at(ao) which 
may be solved iteratively. In ref. |23j the IR renormalon asymptotics with a e ff in ([I]) 



replaced by the only truncated a t C|TT|) has been considered in details. This calculation 
turns out to be surprisingly sophisticated. For any finite number of iterations in the 
transcendental equation (pT|) the asymptotics contains some new functional dependence 
on N compared to the usual IR renormalon (|7|) . The generally accepted result is restored 
only after performing the infinite number of iterations in (|TT|). However, we will skip 
now the discussion of the role of the first few terms of the series (0) and concentrate 
our attention on f3k with very large k. In this case it will be enough to make only one 
iteration in ([□]) (namely ln(at/ao) = — ln(l — b aox) in (|TTD ). 

Consider first the UV renormalon. Moreover, let us restrict ourselves on the UV 
renormalon for pure QED. The first important observation (see [O, 1231) is that the 



series of corrections to renormalon generated by ( j~2]) although have no any reasonable 
small parameter but also is not the asymptotic series. In order to show that it is so 
consider for the moment the simplified version of flTTI , ! 1 "^) with b\ = . Now one has 
immediately for the UV renormalon @ 

1 oo / \ m +l a 

Y—l f3 m i ^ e- y dy = Ya N (-b ) N N\ , , . (14) 

We see that though (3 m themselves are determined by the UV renormalon and therefore 
have the form (3 m ~ m 7 (— bo) m m\ the two most dangerous factors from the (3 m , the ml 
and {—bo)" 1 , have been removed from the series. 

Up to now we have chosen rather arbitrarily the overall normalization of renormalon. 
In order to go further in understanding the role of high order corrections to UV renor- 
malon chain we have to specify the normalization. Consider following ref. [|T^j (and 
almost everybody others in the renormalon business) the correlation function of two 
electromagnetic currents (see also the eqs. (1,2) of ref. JT2J connecting this quantity 

With R e + e -^hadrons) 

U lw =ij dxe^(0\T{j,(x)j u (0)}\0) = (q^ - g, u q 2 MQ 2 ) . (15) 

We will be naturally interested in H(Q 2 ) in the Euclidean domain. 
The simple calculation of the diagram fig. la for k Q gives 



(16) 
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Here the first term in brackets is the parton model prediction, while the integral is ex- 
pected to generate the asymptotic series of the perturbation theory. One may substitute 
into (|16D the one loop running coupling constant and easily reproduce the "traditional" 
(before |12], |22fl ) UV renormalon P5fl . By dots in ( p!6[ ) we have denoted the rest part of 
the perturbative series which is not included into the leading UV renormalon (and is 
expected to be much smaller than the UV renormalon). 

The polarization operator II is trivially connected with the /^-function for the effective 
charge (pf) 



4tt Q z 



(17) 



However, because we are interested only in the asymptotics we may find directly from 
the equation (H) that 



d 



d\n(Q 2 

Thus for our purposes 



-n\(-b a) n = (n + l)\(-b a) 



n+l 



Oi 1 - 



n 



47rU(Q 2 )=J2bnHQ)r 



18) 



(19) 



Now we have to substitute the series for the effective charge a(k) (0) into the integral 
(|16|) and expand the result in series in a(Q). The effective method which allows to find 
the coefficients of such expansion was developed in [23|. First of all it is convenient to 
introduce the new variables 



-b a 



-bi/b 2 



(20) 



P 

For QED b = — , b\ = — ^ and f3 = -^j. With this new variables the truncated 
effective charge in the first nontrivial approximation takes the form 

1 a . ( k 2 ^ 



bo 1 — ax + Pa ln(l — ax) 



x 



In 



,Q 2 



The contribution of the n-th term of the formal expansion (|1 
operator (|Ti]) is of the form 



(21) 

to the polarization 



n, 



12vr 2 (-6 )" 



-i n+l 



ax + /3aln(l 



ax) 



xe 



€ dx 



(22) 



Note that this is the n-th term of the expansion in the series in truncated running 
coupling a t (k) and it still contains the whole series in a(Q). As before we are looking 
for the iV-th term of the series in a(Q). Below we will see that the most important 
contribution comes from n ~ N. The simple binomial expansion now gives 



n, 



-Pn 



127T 2 



a 



1 — ax 



-Pn 



12tt 5 



-b n ) n 



"+ 1 (n + m)\ 

m=0 

n 



—[3a ln(l — ax) 



a 



ax 



i n+l 



E 

m=0 



mini 
n + m)\ 



rami 



1 — ax 
-(3a ln(l 



xe x dx 



ax) 



ax 



e~ x dx. (23) 
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Here in the second equality we have used that 

ax 1 



1 



ax 



ax 



(24) 



and than have neglected the —1 (one may easily see that this less singular contribution 
will lead to the ~ 1/N correction to the asymptotics). 

For estimation of the iV-th order contribution of perturbation theory we will use the 
formula for iV-th term of the expansion of the integral in powers of a 



e- l dt 
1 — at 



(3a 



1 — at 



In 



1 - at 



a N Nl^- 
kl 



(25) 



N 



Here both m and k are supposed to be large. In order to derive (|25|) one has to use the 
asymptotics of gamma-function together with the trivial identity 



Hp)Y 



lim { ^~ 



(26) 



Formula (^) allows now to calculate the iV-th term of the expansion of (|23|) in series 

N 




in an (a N , a = — 6o«o 



{n n } 



-Pn 



N 



127T 2 (-6 )" 



my: 



(3 m (n + m)\ 



127T 2 



-boa ) 



N 



q {n + m)\ mini 



In 



N 



n + m 



(—b ) n nl V n 



N 



(27) 



Here we have neglected m compared to n in the argument of logarithm because effectively 
m ~ ln(iV) while as we will see in the moment n ~ N. Combining together (^,^) and 
the asymptotics of j3k (H) one finds the equation for b N 



37t6 ^ N n {-b ) n nPnl 



(28) 



Here (3 = —b\jb\ = 9/(4JV/) for QED. This equation is even further simplified after 
substitution 



b N 

c N 



-b ) N N?Nlc N , 



(29) 



n <N n(-3nb ) 



N 



E 

/ n<N 



n 



1 + 



The solution to this last equation is evident cn = const x N 1 ^? and for b^ one has 

b N = const{-b ) N N^ + ^J Nl . (30) 

So we arrived at the surprising result: taking into account all possible insertions to 
the renormalon chain allowed to change the power of N in the asymptotics, which for 
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many years was thought to be determined by only the two terms of the Gell-Mann-Low 
(3 -function 6q and b\ . 



The overall constant in fl30|) could not be found in closed form. In terms of the 
equation (|29|) this constant is determined by the initial condition at N ~ 1. But for 



N ~ 1 the equation (|29|) is no more valid as we have indicated explicitly by writing 
(l + 0(l/n)). 

This reduction to the initial value problem shows how one may reformulate the 
problem of the overall normalization of renormalon. Finally, the solution looks almost 
like tautology. One is allowed to look for the renormalon asymptotics in the form 

R N = A n N\-b ) N N\ . (31) 

Here n is the number of terms of the perturbation theory which were calculated explicitly 
(n < N). The normalization constant A n may be found with only the ~ 1/n accuracy. 
This is the important difference between the renormalon and instanton |26| induced 
asymptotics. For instantons not only the overall amplitude of the asymptotics is known 
but also the ~ 1/N corrections to this asymptotics were considered J28| ^] (see also 



However, one may consider the renormalon calculus only as a way to extend the 
explicit perturbative calculations by one more approximate term (and this is the most 
straightforward application of renormalon). In this approach the equation (|3T|) still is 
quite informative. It shows that after the explicit calculation of N terms of perturbation 
theory one will immediately found the N + 1-st term with at least the ~ 1/N 2 accuracy. 



4 UV renormalons with many hard photons 

Our formula (^) for asymptotics of the coefficients b^ would be a nice new result if 
published 3-years ago. However, as we have told in the introduction Vainshtein and 



Zakharov in their preprint of April 1994 [12] have shown that the traditional UV renor- 



malon (|2|) (fig. la) does not give the largest contribution to the asymptotics. They have 
found a series of new diagrams (starting from two three-loop diagrams) which generate 
the asymptotics much larger than (0). The authors of |12[ have used the OPE in order 



to find the contribution to polarization operator of these new diagrams (see Parisi 0, 
who first proposed to use the OPE for renormalon calculus). Finally, the new result for 
n(Q 2 ) reads 

a{k)\ 2+1 Q 2 d*k 



^—Um (32) 

with 7 for the QED case 




2N f 



(33) 
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Now it is easy to substitute the two loop running coupling ( |TT|pTD into (j32|) and find the 
UV renormalon of the form ([31]) with 

*=l + 7 + ^ • (34) 

However, it is clear that this result for 5 will be changed immediately if one substitutes 
the all-loop a e ff into ( |3~^) and repeats the simple calculation which we have performed 
in the previous section for n(Q 2 ) ([l6|) . The magnitude of the correction to 5 (0) now will 
depend on the explicit value of the overall normalization of (p2[) which may be extracted 



from [12, 22 



On the other hand, if we are going to substitute the exact running coupling into 
we should also take into account all possible corrections to this formula 

IW(Q') = cmst l Q {^gf (1 + E ^ M) 9^t (35 ) 

If the coefficients here behave like 

TT m ~ — (36) 

m 

for large m, taking into account of this series will change again the exponent 5 ([34],|3l]). 
Moreover, if n m grows faster with m than fl3T)| ) it will be a catastrophe for asymptotics 
as may be seen from (|29|). The eq. ([36]) is likely the upper bound for the coefficients of 
the expansion 



Unfortunately we do not see now the clear way to estimate the asymptotics of the 
coefficients tt to in fl3"5|). It is rather probable that the corrections to 8 due to the 7r m and 
a e ff{k) 2+1 will compensate each other and the result of refs. |12, ^] will be restored. 



Furthermore, the authors of [T3] have demonstrated in their conclusions that they are 



ready to meet any surprise from the high order corrections to the renormalon. Therefore, 
even if such compensation do not take place and 5 (B3p is changed this result will not 
be in complete disagreement with [|12] . 



Anyway, after one finds the new 5 from fl3"o] ) the overall normalization of the UV 
renormalon will be found with only ~ 1/n accuracy and only after the explicit calculation 
of the first n terms of the series. We write here ~ 1/n although it may be some power 
of it ~ {l/n) k . Moreover, may be just the better understanding of the high order 
contributions to renormalon like in ( [35] ) may help to estimate, with what accuracy one 
may found the n + 1-st term of perturbation theory after explicit calculation of n terms? 



5 High order corrections for the IR renormalon and 
the best accuracy of perturbative QCD 

The IR renormalon have attracted much more interest during the last few years than 
the UV one. This conclusion may be drawn even by simple counting of the number of 
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publications. It seems that most of the authors do not consider the UV renormalon as 
renormalon at all. 



The reason for such asymmetry is quite evident. In QCD the UV renormalon is Borel 
summable while the IR renormalon is not. Physically this non-Borel-summability means 
that depending on the details of chosen resummation procedure the obtained predictions 
for observables will vary by some power corrections. Namely for the case of light fermions 
and inclusive cross sections this correction will be of the form ~ (Aqcd/Q) 4 - There is 
absolutely no way to find the amplitude of this nonperturbative ~ Q~ A contribution 
from the perturbative renormalon calculus. 

If there is no rigourous way to calculate the ~ (Aqcd/Q) 4 correction one may try at 
least to extract them from the comparison of theoretical prediction with the experiment. 
Naively one may take say the experimental value of R e + e -^hadrons subtract the parton- 
model contribution and 2 — 3 known ~ a s corrections and look for the power corrections. 
However, this procedure certainly will not work at least for sufficiently large The 
more or less reasonable procedure is the following: one has to calculate a huge number 

Nm(Q) = — (37) 
b a(Q) 

terms of perturbation theory (Nm is the function of Q !) and subtract them from the 
experimental ratio R(Q). The rest will be the needed (A/Q) 4 correction. 

The first problem (if not worry about the terrible analytical calculations!) which will 
encounter us upon performing this procedure is the UV renormalon. The series for UV 
renormalon blows up at N = l/(b^a(Q)) - twice before the critical value (|37|). However 
the way to avoid this problem seems to be rather clear. As we have seen ([E],|]) the UV 
and IR renormalons originate from the very different regions of variation of the internal 
momentum in the diagram (k in the fig. 1). Naturally one may divide the integral into 
two parts k < Q and k > Q and than obtain the result in the form 

n = J2(-b a) N N\ + £ ( b i*YN\ . (38) 

UV IR,N<N IR \ Z / 

Here the first series is much larger than the second one but allows the explicit (Borel) 
summation. So one has to sum up exactly the series for the UV renormalon in (^) and 
then subtract this resummed contribution from the experimental value of R e +e-->hadrons- 

Our goal is to reach the best accuracy of the perturbative QCD ~ (A/Q) A . There- 
fore the summation of the UV renormalon in (38) also should be done with ~ (A/Q) 4 



accuracy. This is also not so easy to do because even the smallest term of the series 
for the UV renormalon is of the order of ~ (A/Q) 2 - much larger than the accuracy we 
want. For example, if one simply substitutes the 1-loop a(k) into naive renormalon ([16]) 

1 the success of the QCD sum rules approach at low (~ lGeV) energies is usually considered as the 
indication that the actual power corrections in QCD have some additional enhancement compared to 
the typical uncertainty of renormalon resummation. However, this enhancement will became less and 
less important at high energies which only may be a subject of perturbative QCD. 
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the accuracy of the resummed value will be only 0(a 2 ) and so on. It is quite probable 
that in order to sum up the UV renormalon with the (A/Q) 4 accuracy one still has to 
calculate exactly about N IR = 2/(b a) terms of the series. 

However, suppose that accurate enough summation of the UV renormalon in (|38|) has 
been done. Let us see, what can we say about the IR-renormalon part of the polarization 
operator in view of our experience of working with high order corrections to renormalon 



chain? The corresponding contribution to II ( |17|{i9|) now has the form of the integral 



([[]). Simply repeating the calculations (^2|-^) one gets 



N 



But p n ~ b n now is determined by the first UV renormalon and thus (3 n ~ (— bo) n n\ . 
Therefore we see that the sum of the multiloop corrections to the IR renormalon leads 
to the correction of the relative order of magnitude 

~ V - f" , ~ E(~2) n = oo . (40) 

The series of the corrections to renormalon is not asymptotic, but is ugly divergent. It 
looks like some interference of the IR and UV renormalons. 



Nevertheless, as we will see now the divergency itself of the series (f|0]) is not con- 
nected with the infrared physics. As we saw in (^) the N\ in the leading contribution 
to the IR renormalon appears simply due to the N-th power of the large logarithm 
\n(Q 2 / k 2 ) N . Therefore as we have told after the eq. (Q) this trivial renormalon 



x 



N 



contribution falls down into the infrared momentum region like k 2 ~ Q 2 e~~ N l 2 . On the 
other hand, as we have seen in the section 2 for high order corrections to the renormalon 
chain the same iV! appears not only due to the power of logarithm but also due to the 
more and more complicated combinatorics. This change in the origin of the N\ would 
not be taken into account if the series of corrections to renormalon was convergent. 
However, for divergent series of the kind of fl40| ) one should naturally to reformulate 
the method of calculation of IR renormalon in order to be able to control from what 
distances each contribution came. To this end let us rewrite the renormalon ([]]) in terms 
of the series in the powers of \n(Q 2 /k 2 ) 

k 2 dk 2 „ x ^ rk 2 dk 2 +1 




= J2a(Q) n+1 F n (a(Q)) 

n 

Now all the high order corrections to renormalon chain are hidden in F n (a). The function 
F n (a) itself may be expanded in the asymptotic series. But it would be the Borel 
summable series. 

In the previous two sections we have tried to show that there is a lot of open questions 
concerning the UV renormalon. Therefore, now we will illustrate only by the simple toy 
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example how the functions F n (a) ( fill) will be found after one solves all this "ultraviolet 
problems". Suppose the running coupling a(k) in ( fiTf ) is described by the simplest Borel 
integral 



a 



«(*)=/ urn e-«dy = aJ2(-bo&) n n\ , (42) 
Jo 1 + b a{x)y n=0 



Here a(x) is simply the one loop running coupling for the momentum k. The trivial 
expansion in the series over x now gives 

dy (43) 



;i + b a y) 



n+l 



Of course this F n (ao) may be expanded in the asymptotic series in ao- For very large 
number n the integral in (fiSJ) is even further simplified. For example for the last F n -s 



which one is allowed to work with (with n m Njr (|37f) ) the expression (f43f) reduces to 
almost trivial result Fn ir = 1/3. 



6 Event shape renormalons and conclusions 

The IR renormalons for the event shape variables have became the popular topic of 



high energy physics in last couple of years (see e.g [2(| and references therein). It was 
recognized that among the experimental values, which can not be directly related to the 
usual operator product expansion in Euclidean domain, one may find the quantities for 
which the soft gluon contribution is much less suppressed than in ([[]). In particular, the 
corresponding integral may even contain only the first power of k instead of the usual 
fourth power (JJ) 

Res = I <*(k)% . (44) 

It is clear, that such contribution will generate the very large ~ A/Q nonperturbative 
correction. Like for other renormalons this correction currently can not be found within 
ab initio analytic calculation. One may suppose however, (as we have done in fact in 
(|37|)) that the smallest term of the series associated with (|4"4]) may be considered as 
the proper order of magnitude estimate of this ~ A/Q correction. Then in order to be 
able to extract from the experiment the pure nonperturbative part one should calculate 
explicitly the 

= ai 1 (45) 

terms of perturbation theory. Although, formally this number Nes{Q) still is paramet- 
rically very large ~ 1/a , it is 4 times smaller than the corresponding value for the usual 
IR renormalon (|37|) . 
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Moreover, it is usually assumed |Ty|] that the actual ~ A/Q corrections may be 
sufficiently enhanced compared to the smallest term of the series. Suppose that this is 
just the case and the nonperturbative A/Q contribution is in C times enhanced compared 
to the smallest term with some C>1. In this case one may easily modify the eq. ([45]) 

1 



N ES (Q) 



2b a(Q) \ b a(Q) 



^ (46) 



While writing this formula we still have assumed that bo<yln(C) <C 1. Naturally the 
difference ANes between (]46D and (|45[) became relatively less important for smaller 

coupling ANes/Nes ~ Jaln(C). 

Now let us return to the main subject of our paper and consider the high order 
corrections to the Event-Shape renormalon. To obtain the Event-Shape renormalon 
([44]) one has to consider some rather special characteristic of the jet distribution (thrust 
for example in ]17|])- Therefore, it is natural to suppose that the high order corrections 



to the effective coupling in ( [441) are still gathered by the traditional renormalons, and 
especially by the UV one. Under this assumption, repeating again the calculation (^-|28|) 



one gets for the iV-th term of the expansion of ([44]) 



{REs} N ~(2b a) N N x N\\^r m \ , (47) 



,m=l 



where the multiloop corrections r m are expected to behave like 



{2b ) m m\ 



-l/2) m . (48) 



Thus we see that the series of multiloop corrections to the Event-Shape renormalon 
will naturally be described by the simple summable geometrical progression. However, 
this stability of the Event-Shape renormalon also reflects the fact that the best accu- 
racy which in principle may be achieved within the pure perturbative calculation for 



(44) is much worse than for the usual IR renormalon ( ~ A/Q instead of ~ (A/Q) 



1 



respectively ). 

To conclude, the goal of this paper was mainly to outline the problems which should 
be solved upon taking into account the multiloop corrections to the renormalon chain. 
In particular, we have shown that these high order corrections play very different role for 
UV, IR and Event-Shape renormalons. The renormalon contribution to high order terms 
of perturbation theory for Event Shape variables are most hugely divergent compared to 
other renormalons. However, just the Event-Shapes renormalons turns out to be most 
stable under taking into account the multiloop corrections to renormalon chain. On the 
opposite side, the traditional IR renormalon is much weaker than other renormalons, 
but, as we have shown in Section 5, it is most sensitive to the inclusion of complicated 
insertions into the chain (divergent geometrical progression). Nevertheless, one is still 
able to associate with the IR renormalon some perturbation-theory-like series. Only the 
coefficients of this series should themselves be the sums of Borel summable series in av 
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The special importance of the traditional IR renormalon in QCD is due to the fact that it 
gives the principal limitation on the accuracy of pure perturbative calculations. Finally, 
the resummation of high order corrections to any type of renormalons is sufficiently based 
on the use of UV renormalon. In Section 3 we have shown, how taking into account of the 
high order corrections to the single UV renormalon chain may be reduced to solution of 
a simple self-consistent equation. A considerable progress have been achieved last years 
in calculating the UV renormalon asymptotics with many hard photon (gluon) chains 



22| , [30|| . Nevertheless, in Section 4 we speculate about the possibility to change 
the asymptotics again by taking into account the multi-loop insertions into this many - 
chain renormalon. 

The most direct use of renormalons, as well as any other asymptotic estimates, will be 
to predict the value of n + 1-st term of perturbation theory after the explicit calculation 
of n terms. However, contrary to the case of Instantons f26|| , the overall amplitude of 



renormalon itself can not be found without accurate multiloop calculation. Moreover, 
even if the asymptotics of perturbation theory will be found, the corrections to it may 
be also very important. For example in QCD |29| the corrections to the asymptotics 
are of the order of ~ N^/N which makes practically useless the asymptotic estimates in 
this case. Another application of renormalons may be to use some kind of renormalon 
resummation in order to extract the pure nonperturbative-QCD contributions from the 
real high-energy experiment. However, in this case one may hope either on the possible 
(occasional) huge enhancement of the nonperturbative corrections, or to wait a few 



decades before the needed number of terms of perturbation theory ( |37| , ^5| , ^6D will be 
calculated explicitly. 
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